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Abstract: 

The air fluorescence detectors (FDs) of the Pierre Auger Observatory are vital for the determination of 
the air shower energy scale. To compensate for variations in atmospheric conditions that affect the energy 
measurement, the Observatory operates an array of monitoring instruments to record hourly atmospheric 
conditions across the detector site, an area exceeding 3,000 km 2 . This paper presents results from four 
instruments used to characterize the aerosol component of the atmosphere: the Central Laser Facility 
(CLF), which provides the FDs with calibrated laser shots; the scanning backscatter lidars, which operate 
at three FD sites; the Aerosol Phase Function monitors (APFs), which measure the aerosol scattering 
cross section at two FD locations; and the Horizontal Attenuation Monitor (HAM), which measures the 
wavelength dependence of aerosol attenuation. 



Introduction 

The Pierre Auger Observatory in Malargiie, Ar- 
gentina employs four fluorescence detector (FD) 
telescopes to obtain calorimetric estimates of air 
shower energies. The atmosphere, which acts as 
the calorimeter, is constantly in flux, so aerosol 
conditions are measured hourly at each FD loca- 
tion and stored in an offline database for use in the 
reconstruction of showers. 

The aerosol parameters of interest are those that 
affect light attenuation and scattering: a(h), the 
aerosol extinction coefficient, and r(h), the verti- 
cal aerosol optical depth (VAOD), which is the in- 
tegral of a(h) from the ground to altitude h; P(0), 
the normalized differential scattering cross section, 
or phase function, as a function of scattering an- 
gle 8; and the wavelength dependence of aerosol 
scattering. We describe measurements of r{h) and 
its wavelength dependence by the CLF, lidars, and 
HAM, and observations of P(6) by the APFs. 



Optical Depth Measurements 
Central Laser Facility 

The CLF produces calibrated laser "test beams" 
from its location in the center of the Auger surface 
detector [4, 9]. Located between 26 and 39 km 
from the FDs, the CLF uses a pulsed laser oper- 
ating at 355 nm. The atmosphere scatters nearly 
equal amounts of laser light toward each FD eye, 
where it is recorded by the FDs. With a nomi- 
nal energy of 7 mJ per pulse, the light produced 
is roughly equal to the amount of scintillation light 
generated by a 10 20 eV shower. 

Among the many measurements provided by the 
CLF test beams are hourly observations of a(h) 
and t(/i). The laser fires sets of 50 shots every 
quarter-hour, and the track profiles are averaged to 
obtain hourly aerosol measurements. The CLF es- 
timates the aerosol content using an iterative pro- 
cedure that does not require absolute photometric 
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Figure 1: CLF VAOD measurements at Los Leones Figure 2: VAOD at 4.5 km. Los Leones, at an alti- 
and Coihueco, 2004 - 2006, showing mean VAOD tude ~ 280 m below Coihueco, observes a systemat- 
profiles and their 68% confidence limits. ically larger VAOD. 



calibrations of the FD or laser. The procedure be- 
gins by normalizing hourly average laser profiles 
L(h) by an "aerosol-free" reference profile L ie f(h) 
measured during extremely clear conditions. The 
initial VAOD estimate is then 

= _ lnL(h)-lnWfe) 
1 + csce(ft) 

where e is the elevation angle to the track point at 
altitude h. An estimate on (h) comes from the slope 
of the VAOD, and is used to iteratively correct the 
light profile for aerosol scattering. At the end of 
the calculation, the final a(h) is normalized by the 
VAOD at a height h c above the bulk of aerosols: 

j " a(h)dh = Ti(h c ) (2) 

" "■ground 

The final VAOD is found by integrating a(h). Fig- 
ures 1 and 2 depict the VAOD distribution recorded 
at Los Leones and Coihueco between 2004 and 
2006. Typically, the VAOD at 4.5 km is 0.03, with 
statistical uncertainties of ~ 0.01. 

The uncertainties in each VAOD measurement are 
dominated by statistical fluctuations in the hourly 
average light profiles, although systematic effects 
due to the FD and laser relative calibrations and 
the choice of aerosol-free reference nights are also 
present. The last effect has been checked with a 
separate analysis method that uses CLF laser sim- 
ulations and a simple two-parameter exponential 



model of the aerosol density. The results, which 
are independent of clear-night calibrations, closely 
match the standard CLF estimates (Figure 3). 

The CLF can also detect clouds, which appear as 
sharp steps in the VAOD profile. When a cloud is 
present, the lowest base height ft-basc is recorded. 
For heights h < hbase> an d T {h) aie consid- 
ered valid for air shower analysis. 

Elastic Backscatter Lidars 

In addition to the CLF, the observatory employs 
three scanning elastic lidar stations, with a fourth 
under construction, to record a(h) and r(h) at 
every FD site [2]. Each station has a steerable 
frame mounted with a pulsed 351 nm laser, three 
parabolic mirrors, and photomultiplier recorders. 
The station at Los Leones also includes a separate, 
vertically-staring Raman lidar test system, which 
can detect aerosols and the relative concentration 
of N2 and O2 in the atmosphere. Since the lidar 
hardware and measurement techniques are inde- 
pendent of the CLF, the two systems have essen- 
tially uncorrected systematic uncertainties. 

Every hour, each lidar sweeps the sky in a set pat- 
tern, pulsing the laser at 333 Hz and observing the 
backscattered light with its optical receivers. The 
sweeps occur outside the FD field of view to avoid 
triggering the detector. By observing the backscat- 
tered light, the lidar can determine a(h) and r(h) 
using a straightforward numerical inversion [5], As 
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Figure 3: VAOD at 4.5 km from the CLF in different Figure 4: VAOD at 4.5 km at Coihueco, Oct 2006 
calibration epochs, compared to laser simulations. Jan 2007, observed by the CLF and Lidar. 



shown in Figure 4, the lidar and CLF results are in 
reasonable agreement, despite large differences in 
operation, analysis, and viewing regions. 

The lidar is also well suited to detect cloud layers, 
which create significant echoes in the backscat- 
tered light signal; these can be automatically de- 
tected by a simple gradient/threshold algorithm. 
Under certain circumstances, the lidar can also es- 
timate cloud optical depths [2]. The lidar is cur- 
rently accumulating an extensive hourly database 
of cloud height, sky coverage, and optical depth. 



/ and g determine P(9), and are affected by local 
aerosol characteristics. 

In the Auger Observatory, / and g are measured 
by Aerosol Phase Function monitors (APFs) lo- 
cated several km from the FDs at Coihueco and 
Los Morados [3], Using a collimated xenon flash 
lamp, each APF fires an hourly sequence of 350 nm 
shots horizontally across the FD field of view, cov- 
ering 30° to 150° in azimuth. The scattering pa- 
rameters / and g can be determined simply by fit- 
ting the horizontal light track recorded by the FD. 



Scattering Measurements 

The FD reconstruction must not only correct for 
the attenuation of scintillation light, but also sub- 
tract Cherenkov light scattered into the FD field of 
view. Therefore, the scattering properties of the 
atmosphere should be well-understood. Aerosol 
scattering is highly nontrivial and depends on the 
physical properties of the aerosols, but its distri- 
bution as a function of scattering angle 8 can be 
approximated by the simple parameterization 



P(9) 



1-g 2 

47T 



(3) 



3 cos 2 e - 1 



(l+.g 2 -2 ff cos6») 3 / 2 J 2(l + 5 2 ) 3 / 2 



where g = (cos 9) measures the asymmetry of 
scattering and / determines the relative strength of 
forward and backward scattering. The quantities 
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Figure 5: Scattering parameter g measured 
Coihueco between June 2006 and March 2007. 
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Ten months of APF measurements at Coihueco 
have yielded a site average of g — 0.59 ± 0.08 
for the local asymmetry parameter, excluding clear 
nights when g = 0. The distribution of g, shown in 
Figure 5, is comparable to measurements reported 
in the literature for similar climates [1]. 
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Wavelength Dependence 

The attenuation of light by aerosols is expected to 
have some wavelength dependence, and this is typ- 
ically parameterized by a power law 

t(A)=t„.(V) 7 (4) 

In this expression, To is the optical depth measured 
at the reference wavelength Ao = 355 nm, and 7 
is the so-called Angstrom exponent of the depen- 
dence (7 « 4 for molecular scattering). 
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Figure 6: HAM fit to CCD response after molecu- 
lar subtraction for one measurement set. 
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Figure 7: Distribution of the Angstrom exponent 
observed by the HAM, July 2006 - February 2007. 

At the Auger Observatory, measurements of 7 
are performed by two instruments: the robotic 
telescope FRAM, described in detail in [8]; and 
the Horizontal Attenuation Monitor (HAM). The 
HAM consists of a high intensity discharge lamp 
located at Coihueco, viewed by a CCD camera 
placed ~ 45 km away at Los Leones. Using a filter 



wheel, the camera records the aerosol extinction 
coefficient, and hence the aerosol optical depth be- 
tween the two sites, at five different wavelengths. 

Figure 6 shows a typical HAM fit used to estimate 
7. The uncertainties are dominated by measure- 
ment fluctuations, and include a systematic effect 
due to subtraction of the estimated molecular opti- 
cal depth between Los Leones and Coihueco. The 
distribution of observed 7 values, plotted in Fig. 7, 
is consistent with other measurements and physical 
expectations [8, 6]. 

Impact of Aerosols on Shower Measure- 
ments 

The measurement uncertainties in the VAOD, 
phase function, and wavelength dependence of 
aerosol attenuation have been propagated in the 
reconstruction of a select number of high-quality 
air showers [7]. The dominant uncertainty comes 
from the VAOD, which contributes 5.5% to the un- 
certainty in shower energy, while the phase func- 
tion and wavelength dependence both contribute ~ 
1% to the uncertainty in the energy. The effect of 
VAOD observations on X max , the depth of shower 
maximum, is 4 g cm~ 2 , while the wavelength de- 
pendence and the phase function contribute uncer- 
tainties of 1 and 2 g cm~ 2 , respectively. 
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